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Abstract 

The wormlike chain model of stiff polymers is a nonlinear a-model in one spacetime 
dimension in which the ends are fluctuating freely. This causes important differences 
with respect to the presently available theory which exists only for periodic and Dirichlet 
boundary conditions. We modify this theory appropriately and show how to perform a 
systematic large-stiffness expansions for all physically interesting quantities in powers of 
L/£, where L is the length and £ the persistence length of the polymer. This requires 
special procedures for regularizing highly divergent Feynman integrals which we have de- 
veloped in previous work. We show that by adding to the unperturbed action a correction 
term .A corr , we can calculate all Feynman diagrams with Green functions satisfying Neu- 
mann boundary conditions. Our expansions yield, order by order, properly normalized 
end-to-end distribution function in arbitrary dimensions d, its even and odd moments, 
and the two-point correlation function. 



1 Introduction 

Recently, the study of biopolymers has become a subject of increasing interest in the research 
of biological materials. Forming networks, constituent filaments play an important role in 
the structure and function of living cells and other biological entities. 1 ' 2 Recent advances in 
visualizing and manipulating single cytoskeletal filaments 3 ' 4 and DNA 5 ' 6 have inspired the study 
of conformational characteristics of biopolymers in single- molecule experiments. 7 ' 8 

Biopolymers may be flexible, such as DNA, semiflexible or stiff, such as actin, or rigid, 
such as microtubuli. For a polymer with stiffness k, the persistence length of tangent-tangent 
fluctuations is £ = 2k/ {d — 1)/cbT. Typical values of £ in biopolymers range from several nm 
to a few mm. Some example are DNA, £ « 50 nm 5 ; spectrin, £ « 10 nm 4 ; actin, £ « 17/mi, 
and microtubules, £ ~ 5 mm. 3 

The length L of a polymer is conveniently measured in units of £, which defines a dimen- 
sionless parameter characterizing the inverse stiffness, or the reduced length I = L/£, which we 
shall also call flexibility. For high flexibility I, a polymer behaves approximately like an ideal 
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random chain with freely joint links. 9 For small flexibility /, is behaves like a rigid rod. Un- 
der the influence of longitudinal compressing forces these filaments exhibit the classical Euler 
buckling. 10 

Most constituent filamentary biopolymers are semiflexible with intermediate values of Z. If Z 
is of the order of unity, the stiffness is important but not overwhelming. Effects of stiffness have 
been observed even for large flexibilities up to I m 1000, for instance in long strands of DNA. 
Stretching these at large forces has shown significant deviations from a random chain model. 6 
It is therefore important to develop a reliable theory to calculate the influence of stiffness upon 
polymers for the entire range of flexibilities. 1,2 

If a polymer has a sizable stiffness, self-avoidance effects can be neglected due to the energetic 
suppression of configurations where the filaments fold back onto themselves. An appropriate 
model for a theoretical description of semiflexible polymers is the wormlike chain model pro- 
posed by Kratky and Porod in 1949. 11 This model is formulated most directly in terms of 
path integrals. The calculation of the statistical properties has mostly been based so far on 
the equivalent diffusion equation. 12 For a recent review see Yamakawa's textbook, 13 the review 
article by Chirikjian and Wang, 14 and Chapter 15 of the textbook of one of the authors. 15 

The wormlike chain model is the minimal model of a polymer with arbitrary stiffness. A 
central feature of this model is the local inextensibility of filaments which is mathematically 
accounted for by constraining the length of all tangent vectors to unity. Being a sequence of 
o?-component unit vectors, a polymer is equivalent to a particle moving on the surface of a unit 
sphere in d dimensions, which is described by a quantum- mechanical 0(oQ-symmetric nonlinear 
(T-model in one spacetime dimension. 

Another equivalence exists to the continuum limit of the classical Heisenberg chain of unit 
magnets, In fact, the partition function of the lattice model was calculated 35 years ago by 
Stanley 16 in d = 3 dimensions. In the magnetic context, however, the quantities of interest in 
polymer physics were never considered. 

Because of the nonlinearity of the wormlike chain model, only a few of the statistical prop- 
erties of the stiff polymer can be calculated analytically, the most prominent being the mean 
square ( R 2 ) of the end-to-end distance R, and a number of higher moments (R 2n ), for n < 14. 
The lowest even moments have been calculated a long time ago analytically 17 and numerically 18 
by solving recursively the diffusion equation on the unit sphere in three dimensions. The cal- 
culations have recently been extended for d = 3 up to n = 13 in Ref. 19 and for arbitrary d 
up to n = 14 in Ref.. 20 By this technique it is impossible, however, to derive the end-to-end 
distribution function for all persistence lengths. To derive this central quantity of the wormlike 
chain, one has to resort to perturbation schemes. So far, only approximate but highly accurate 
expressions for the end-to-end distribution of two- and three-dimensional stiff polymers were 
obtained for all persistence lengths in Refs.. 19-21 

The end-to-end distribution function P(R; L) is a principal physical quantity characterizing 
the statistical properties of a single polymer of length L. For models like the wormlike chain 
which has only short-range interactions between monomers, it gives also the probability density 
of finding any two monomers at spatial distance R = x(s) — x(s'), not just the end points, if 
we identify L with the distance \s — s'\ between monomers measured along the chain. 

In the random-chain limit of large /, the distribution P(R; L) is known exactly 13 ' 15 and 
can be approximated by a simple Gaussian. For moderately large /, the corrections to the 



2 



Gaussian distribution were calculated up to the order I 2 in three dimensions by Daniels, 22 and 
for arbitrary d by one of the authors. 15 

For polymers close to the rod limit, all moments and the distribution function P(R; L) 
were calculated as expansions in powers of / for d = 3 in Refs.. 17 ' 23 ' 24 In the same limit, the 
distribution function has recently been found as an infinite series of parabolic cylinder functions 
for d = 2, and of Hermite polynomials for d — 3 in Ref.. 25 The generalization to d dimensions 
is given in Chapter 15 of the textbook. 15 

These distribution were derived from the path integral of a simple harmonic oscillator. 25 For 
large stiffness, they are in good quantitative agreement with Monte Carlo data. They disagree, 
however, with the exact expansions of Ref., 23 except for a few lowest orders. Since the failure 
has a perturbative character, it cannot be compensated by a proper normalization. The reason 
is that the harmonic path integral approximates correctly the wormlike chain model only at 
lowest flexibility. It certainly needs higher-order corrections. 

If we want to calculate such corrections in the path integral formalism, the choice of phys- 
ically appropriate boundary conditions is essential. The open-end boundary conditions of the 
wormlike chain model are quite tedious for the perturbative calculation of path integrals. For 
a harmonic oscillator, they can be replaced by the simpler Neumann boundary conditions as 
proposed in Ref.. 25 However, to describe the properties of a wormlike chain, the harmonic fluc- 
tuations must be suppressed at the endpoints. This can be done by incorporating an important 
correction factor into the path integral with Neumann boundary conditions. 

Calculation schemes based on periodic boundary conditions were considered in Refs.. 26 ' 27 
Path integrals for modified wormlike chain models with softened inextensibility condition have 
been suggested in Refs.. 28-31 These models, however, can describe the behavior of the stiff 
polymer only roughly. Some conformational properties were calculated in Refs.. 32,33 Addi- 
tional insights into the properties of the end-to-end distribution function of polymers with 
intermediate stiffness have also been provided recently by Refs.. 34 ' 35 

Despite a number of recent developments, the theoretical understanding of the statistical 
properties of a semiflexible chain in isolation remains much less developed than for the flexible 
chain, and experimental data are often interpreted within the theoretical framework established 
for random chains or rigid rods, respectively. It is therefore important to find an unified analytic 
approximation procedure which will yield a reliable end-to-end distribution over the entire range 
of stiffnesses. In particular, the crossover region between the random chain properties at low 
stiffness and the rigid-rod limit should be described properly. 

In this paper, we develop a systematic path integral approach to find the properties of 
the wormlike chain model near the rod limit. We shall set up a perturbation theory for a 
large-stiffness expansion in powers of the flexibility I and show how to calculate all expansion 
coefficients analytically. All moments, the radial end-to-end distribution function, and the two- 
point correlation function will be obtained in arbitrary dimensions d. By this unified approach 
we reproduce correctly the results of Ref., 23 and show how to go systematically beyond the 
approximate large-stiffness result of Ref. 25 for a polymer in d = 3 dimensions. 
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2 Wormlike Chain Model 



In the wormlike chain model, the polymer is described by a smooth curve x(s) in rf-dimensional 
flat space, with the components x l (s), i = 1, . . . , d. The parameter s G (0, L) is the arc length 



along the curve defined by ds = y(rfx) 2 . In this parametrization, the derivatives dx(s)/ds = 
u(s) are tangent vectors of unit length: 

|u(*)|=l. (1) 

This important property accounts for the local inextensibility of a polymer. All vectors u(s) 
lie on a surface of a unit sphere in d dimensions which has an area Sd = 2ir d / 2 /T(d/2) and a 
constant scalar curvature R — (d — l)(d — 2). 

The local curvature of the polymer is k(s) =\ u(s) |, where u(s) = du(s)/ds. The bending 
energy of the wormlike chain model is quadratic in k(s): 



where we have introduced the inverse stiffness e = k B T j k in units of k B T. The flexibility I is 
simply related to e by I = eL(d — l)/2. 

The statistical properties of the model are completely determined by the partition function 
Z = S, 1 J d d u b 5(u 2 b -I) J d d u a 5(u 2 a - l)z(u b ,u a ) , (3) 
where z(u b , u a ) is a partition function density defined by the path integral 



/■u(L)=U(, / 1 rL 

z(u b ,u a ) = e eLR ' & / V d u(s)5[u 2 (s) - 1] exp -- / dsii 

ju(o)=u„ \ Ze Jo 




(4) 



This coincides with the path integral for the imaginary-time evolution amplitude (ubL|u a 0) of 
a particle on the surface of the unit sphere in d dimensions which possesses an obvious 0(d)- 
symmetry. The length parameter s plays the role of the imaginary time or pseudo-time, and 
the bending energy (2) in the exponent of (4) is proportional to the euclidean action A of a 
particle: E c [u] = k B TA[u]. We shall often use this analogy in the sequel. 

The delta-functional in the integrand of (4) accounts for the local inextensibility (1). At 
the endpoints s — 0, L, this constraint is also accounted for by the two ordinary (^-functions 
in Eq. (3). The restricted ordinary integrals over the initial and final u- values lead to the 
partition function with open ends (3). The path integral (4) for the density is calculated 
with fixed ends, i.e., with Dirichlet boundary conditions (DBC). The extra covariant prefactor 
e eLR/8 j g necessar y to reproduce the correct diffusion equation on the unit sphere 15 ' 38 whose 
Hamiltonian is a Laplace-Beltrami operator. 

The principal physical quantity of a semiflexible polymer is the distribution function 

P(R;L) = /^5 {d) (r- £ dsu(s)j^ , (5) 
with R = x(L) — x(0) being the end-to-end vector. It possesses the moments 



/ / dsds' u(s) ■ u(s') 
Jo Jo 



n/2 



(6) 
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which have been calculated exactly for even n up to high n by solving the diffusion equation 
on a sphere. From this solution we also know the rotationally invariant two-point correlation 
function 

G(s,s') = (u(s).u(s')) (7) 

as being simply 

G(s,s') = eM-\s-s > \/0 (8) 

displaying clearly the exponential decay of tangent-tangent correlations over the persistence 
length £. We shall All expectation values (5)-(7) must be calculable from the path integral 
representation (3) of the partition function by inserting the quantities inside the Dirac brackets 
into the integrand of Eq. (4). 

The inextensibility condition (1) makes the path integral non-Gaussian. For an analytic 
treatment, we may follow two methods. One may either solve the inextensibility condition (1) 
explicitly. Then the 0(<i)-symmetry is realized nonlinearly. In geometric language, 0(d) is 
the isometry of the metric on a surface of the sphere. This method will be used in this paper 
to develop the large-stiffness perturbation expansion of path integrals, in which all physical 
quantities are expressed as power series in /. A second method may be based on enforcing the 
inextensibility condition (1) with the help of a Lagrange multiplier. This leads to calculations 
of the physical quantities as power series in 1/d, to be presented in a forthcoming paper. 40 Both 
options lead to systematic approximation schemes. 

3 Large-Stiffness Expansion 

The large-stiffness expansion of the path integral in powers of e oc I is a special case of weak- 
coupling perturbation expansions, which have recently been developed by us for quantum- 
mechanical path integrals of nonlinear a- models with various boundary conditions. 15 ' 36 ~ 39 

3.1 Calculation with Dirichlet Boundary Conditions 

To make contact with the methods developed in the previous papers, we eliminate the 5- 
functional in (4) by setting u(s) = (q°(s), ? M (s), ) with /j, — 1, . . . , d — 1, and integrating over 
q°(s), which yields q°(s) = a(s) = \Jl — q 2 (s). The partition function density becomes 

z(Qb, Qa) = J dbc V d - 1 q(s) yjg(q(s)) exp j-i- £ dsg^(q) <f (s)<f (s) + eL | j , (9) 

with the metric g^u(q) = 5^ + (1 — q 2 )' 1 q^q u and the determinant g(q) = (1 — q 2 )^ 1 . 
The square root in the invariant measure can be rewritten formally as 

Y[\/g((l(s)) = ex P j^E lo g#(<?( s ))j = ex p|^ ds5(s,s)\ogg(q(s))^ , (10) 
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where the quantity Jq ds5(s, s) = L5(0) counts the infinite number of eigenvalues of the oper- 
ator —5^ u (d/ds) 2 in the space of functions q^(s). Including the exponent of (10) into the action 
we are left with the partition function density 

z(q b ,q a )= I V d - 1 q(s)exp{-A[q]}, (11) 

JDBC 

with the euclidean action 

A[q] = ^ J q L ds [g, u (q] q»(s)q»(s) - e5(s, s) hgg(q(s))} - sL | . (12) 

Repeated indices are summed as usual, i.e., q^q 11 = q 2 . In a perturbation expansion of the path 
integral (9), the inverse stiffness e oc I plays the role of a coupling constant. The large- stiffness 
or small-flexibility expansion is a weak-coupling expansion. The subscript DBC in Eq. (9) 
emphasizes the Dirichlet boundary conditions q^(L) = q£, g M (0) = q%. 

To derive the perturbation expansion of Refs., 15 ' 36-39 it is convenient to rescale the coor- 
dinates q ll (s) — > \feq^(s). This removes e from the quadratic part of the action. The path 
integral for e = is therefore Gaussian and determines the basic free correlation function or 
propagator 

(q»(s)q»(s')) = 5^A(s,s>), (13) 

where A(s, s') is the Green function of the operator — d 2 in the space of functions q(s) with 
appropriate boundary conditions. 

An expansion of the action in powers of e organizes the infinitely many interaction terms. 
They give rise to fluctuation corrections which are calculated by performing all possible Wick 
contractions. These consist of multiple integrals over the pseudo-time s containing products 
of propagators (13) and their s-derivatives. These integrals are mathematically undefined, 
since derivatives of the propagator (13) contain generalized functions or distributions, such as 
Heaviside functions Q(s — s') and Dirac (^-functions 5(s — s'). The problem of defining these 
integrals has, fortunately, been solved in our previous work. 36-38 We have shown that all these 
highly singular integrals are defined uniquely on the basis of the simple physical requirement 
that the path integrals should be independent of the coordinates used to parametrize the unit 
sphere. The new integration rules are valid for path integrals with all boundary conditions. 
Moreover, they are in complete agreement with much more cumbersome calculations in D = 1— e 
dimensions, in which the limit e — > is taken at the end (dimensional regularization) . For one- 
dimensional cr-models with Dirichlet and periodic boundary conditions, our rules have led to a 
perfect cancellation of all UV-divergent diagrams involving powers of 5(0) in each perturbative 
order. 37-39 This permitted us to calculate finite fluctuation corrections to any desired order in 
perturbation theory. 

The method developed in Refs. 15,36-39 can be used now to derive the partition function with 
open ends (3). In coordinates q' 1 (s), this quantity takes the form 

Z = S, 1 J d d -\ b jgj^) J d d -\ a jgjq^)z(q b ,q a ). (14) 

This way of writing Z suggests that we must compute the path integral (9) with Dirichlet 
boundary conditions, and subsequently perform two ordinary integrals over initial and final 
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g M - values. The advantage of this procedure would be that all necessary tools are available from 
our previous work in Ref.. 15 ' 38 We would expand the action (12) in powers of the coordinates 
q^(s) around the straight-line solution connecting the end-points, and evaluate higher-order 
fluctuation corrections with the help of the basic propagator (13) with Dirichlet boundary 
conditions 

A D (s, s') = - | s - s' | /2 + (s + s')/2 - ss'/L . (15) 

The resulting expansion of z(qb, q a ) involves the geometric invariants formed from the curvature 
tensor on g-space. For the unit sphere in d dimensions, only powers of the curvature scalar 
R — (d — l)(d — 2) remain, and the expansion has the form 

e -^ cl [Ag;e] oo 

z(gb,Qa) = (q b L\qafy = — d L k a k (Aq; s) , (16) 

where the classical action in the exponent depends on Ag M = — q a )^ as 

^, [A , ;£] s (M + £ KM!!! + ..., (17) 

and the lowest coefficients dk(Aq;e) have the small-e expansions 

- l + ^(A,) 2 + ^(^H^) [( a,)T + ..., 
0l(A9;E) s E i±zm^A + £ , (^K^_i!^ (Agr + ..., 

. . . Jd-l)(d-2)(5d 2 -17d+18) , s 

a 2 (Aq;e) = e 2K - ^ ^ >- + .... (18) 

From this we can immediately find the partition function (14) with open ends. Inserting 
(16) into (14) yields 

Z = 1, (19) 

to all orders in e. 



3.2 Calculation with Neumann Boundary Conditions 

For a calculation of the important polymer properties (5)-(7) this procedure, although well 
defined, would be too tedious. To save efforts, we prefer to do the calculations from a path 
integral in which the ordinary integrals over endpoints in Eq. (3) are included in the path 
integral by suitable boundary conditions. For a simple harmonic oscillator, these were shown 
in 15 to be the Neumann boundary conditions (NBC), as anticipated by Ref.. 25 These conditions 
may, however, not be used directly for calculating the fluctuation corrections since they are 
physically incorrect. The Neumann boundary conditions admit only fluctuations in which the 
endpoints have zero derivatives g M (0) = q^{L) = 0, thus neglecting all paths with nonzero 
derivatives. We shall see that this omission can be compensated by adding to the action an 
important correction term, so that Neumann boundary conditions may be used after derive the 
correct perturbation expansions. 
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We shall evaluate represent the quantities (5)-(7) from path integrals with Neumann bound- 
ary conditions, by replacing the combined measure of Eqs. (3) and (9) as follows 



[d d -\^lb) fd^qa^a) I V d -\{s) -> / V^qis) J[q b , q a ] , 

J v J v JDBC JNBC 



(20) 



where the Jacobian factor J[qb, q a ] corrects for the missing endpoint fluctuations when restrict- 
ing the paths to zero derivatives at the ends. This factor may be attributed to an extra action 
at the end points 

A co *{q b , Qa] = - log J[q b , q a ] = -[q 2 (0) + <? 2 (L)]/4, (21) 

to be added to (12). 

Our correction term (21) should no be confused with a similar-looking but completely un- 
related term A{[u 2 (0) — 1] + [u 2 (L) — 1]} in harmonic models of stiff polymers, 29-31 where it 
serves to enforce approximately the inextensibility (1) at the endpoints, which in our approach 
is is satisfied exactly for the entire polymer. 

Thus we shall calculate the polymer properties (5)-(7) by performing the averages inside 
with respect to the path integral representation for Z with Neumann boundary conditions 



V d - 1 q(s)exp{-A{q]-A™{q l ,q a ]}. 

NBC 



(22) 



The perturbation expansion of this expression is not straightforward since with Neumann 
boundary conditions, the fluctuations of the free part of the action contain d — 1 zero modes. 
These are due to the d — 1 constant solutions q^(s) = & of the "equation of motion" 
— (d/ds) 2 q^(s) = 0. Their symmetry origin lies in the d — 1 isometries of a sphere in d di- 
mensions. They prevent us from inverting the operator —(d/ds) 2 to obtain a perturbative 
propagator (13). The zero modes must first be extracted from the path integral. 

The proper way of doing this has recently been exhibited for one-dimensional path integrals 
of a nonlinear a-model with periodic paths in Ref.. 39 The method can easily be modified to 
apply to Neumann boundary conditions. To avoid overcounting of the constant paths in the 
elimination process, we have to fix a coordinate system g M (s), for example, by placing the 
center of mass of a polymer at the origin. Let us assume that the end-to-end distance vector 
R = J L dsu(s) points towards the dth direction. In coordinates g M (s), this yields the zero 
average 

L- 1 f L ds<f{s) = 0, n = l,...,d-l. (23) 
Jo 

We therefore introduce the measure of fluctuations without the dangerous d—1 zero modes as 



V' d - 1 q{s)=V d - 1 q{s)5 {d ~ 1) 



dsq^s) 



(24) 



to be used instead of the measure in Eq. (20). 

This can be done most directly for the distribution function (5). Introducing the reduced 
end-to-end distance r = |R|/L, we represent this quantity via the path integral 

P(r;L) = S, 1 f V' d ~ 1 q(s)5 (r - L' 1 ^ ds y/l - q^s)) exp{-A[q] - A™[q h , q a }} . (25) 

JNBC \ JO J 
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In turn, the partition function (22) can be obtained from the distribution (25) as 

roo 

Z = S d drr d - 1 P(r;L) , (26) 
Jo 

where the prefactor Sd reflects the rotation invariance of P(r;L). Since the zero modes are 
already excluded from the path integration in Eq. (25), we integrate simply over r in Eq. (26) 
with help of the 5-function in Eq. (25) to define the path integral (22). This produces a further 
contribution to the action (12): 



A ^ [q] = _( d _ 1) l 0g (L-l ^ ds , 



(27) 



to be referred to as Faddeev- Popov action. For periodic paths this action was found in our pa- 
per. 39 It is a logarithmic Jacobian which compensates the distorting effect of the condition (23) 
on the measure of path integration. This procedure of extracting the zero modes guarantees 
therefore the independence of path integrals of the choice of the coordinate system (23) used 
in the perturbative calculation. 

Thus we rewrite the partition function (22) as a path integral 

Z = J^V'^q^expl-AM-Ariq^qaj-A^iq}} , (28) 

where the paths not only satisfy the Neumann boundary conditions, but have also no zero 
modes. Note that the prefactor Sd in (22), the volume of the isometries, is now absent. 

The expectation values in Eqs. (5)-(7) can now be defined with respect to this partition 
function as follows 

(...) = / NBC ^' d "Vs) (•••) e W {-A[q]-A cor [q b ,q a }-A FP [q}} . (29) 

Note that there is no need to normalize the path integral (29) by Z = Z , since the partition 
function (28) is equal unity as in Eq. (19). This will be verified order by order in perturbation 
expansion of the path integral (28). 

With the definition (29), the radial distribution function is given by the path integral 
P rad (r; L) = (s(r~ L' 1 ^ ds yjl - g 2 (s)^) ^ 

= J NB V >d ^q(s)5 (r - L- 1 ^ ds y/l - q^exp {-A[q] -A™[q b , q a ] -A FP [q]} 
= Sdr^P^L), (30) 

where the Faddeev-Popov action (27) is illuminated by the 5-function, whereas the moments (6) 
of the reduced distance r are found from 

((R 2 /L 2 ) n/2 ) = {r n ) = d 3 y/l - q^J^ 

L- 1 ^ ds^l-q%s)\ e W {-A[q]-A co <[q b ,q a ]-A FP [q]} 
= I V'^qis) exp {-A[q] - A? fa, q a ] ~ A F n P [q]\ = Z n . (31) 
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This looks similar to Eq. (28) except for the last action term which now is 



A F n P [q] = -(n + d-l) log (L- 1 £ dsy/l-q*^ , 



(32) 



rather than (27). The path integral (31) abbreviated by Z n comprises for n = the unit 
partition function (28): Z — Z — 1. 

With the zero modes subtracted, the Green function with Neumann boundary conditions 
representing the lines in the Feynman diagrams acquires the form 



A' N (s, s') = L/3- | s - s' | /2 - (s + s')/2 + (s 2 + s' 2 )/2L , 



(33) 



[ L dsA' N (s,s') = 0. 
Jo 



(34) 



In the following we shall simply write A(s, s') for A' N (s, s'), for brevity. 



4 Partition Function and All Moments up to Four Loops 

We present now the explicit perturbative calculation of the partition function (28) and all 
moments (31) to order e 2 oc I 2 . This requires evaluating Feynman diagrams up to four loops. 
As mentioned before, the fluctuation corrections involve integrals over products of distributions, 
which will be calculated unambiguously using the rules derived in Refs.. 37 ' 38 

For perturbation calculation, we rescale the coordinates g M (s) — > ^/eq fJ '(s) and rewrite the 
path integral (31) as 



NBC 



(35) 



where the total action A^lq', s] = A[q; e] + A cor [qb, q a ', s] + A „ [q; e] reads explicitly, 



A™[q;e] = / ds 
Jo 



a n \og 



l r L 



T^)+\mio g (l-eq 2 ) 



LJo 



dsJl — eq 1 



q 2 (0)+q 2 (L)]-eL 



R 



(36) 



The constant a n is an abbreviation for a n = n + (d — 1). For n = the path integral (35) must 
yield the normalized partition function Z — Zq — 1. 

At e — 0, the total action (36) contains only the kinetic term 

A°[q]= 1 -J o L dsq 2 (s). (37) 
With this free action, the path integral (35) becomes Gaussian which will be defined as 

Z°= [ V' d ~\{s) e- A °M = ( V'^qis) e -W)J L *>?M = l. ( 38 ) 

JNBC JNBC 
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To find the corrections, we split 



(39) 



where the interaction action ^4™%; e] has to be expanded in powers of small e. Its large- stiffness 
expansion starts out as 



A%[q;e]=eA% 1 [q]+6?A™ [?] + ■■■■ 
The first expansion term is 

1 r L 



A^[q] 



2 Jo 



ds {[q(s)q(s)} 2 - p n (s) q~(s 



(40) 



(41) 



where p n (s) = 5 n + [S(s) + S(s — L)] /2 with 5 n = 5(0) — cr n /L. The second expansion term is 
given by 



A^[q] = \[ds[[q(s)q(s)] 2 -\ 



5(0) - 



2L 



\s)) q 2 (s) + ^ ^dsl'ds'qWis') . (42) 



The perturbation expansion of the path integral (35) in powers of e is an expansion in 
terms of expectation values to be calculated with respect to the Gaussian integral (38). For an 
arbitrary functional F[q] of q(s), these will be denoted by 



(F[q]) = [ V'^q{8)F[q\e-W>!>* 

JNBC 



With this notation, the perturbative expansion of (35) reads 



(43) 



Z n = l-(A^[q;e])o + -(A^[q;e} 2 ) -... 



1 



1 - e (A^[q}) + e 2 -(A^[q}) + -«« [?] ) 



(44) 



The expectation values (. . .)o are evaluated by performing all possible Wick contractions with 
the basic propagator (13) and the Green function (33) of the unperturbed action (37). The 
relevant loop integrals 7j and Hi are calculated with our regularization rules when necessary 
and are listed in Appendices A and B. 

We now state the results for various terms appearing in Eq. (44): 



(AT[<A)o 



(d-i) 



^/i + d/ 2 -^A(0,0)-iA(L,L) 



— L— — L {d - l)n , (45) 
8 12 ' v ; 



(A"f[q])o = 



(d 2 - 1) 



2L 



6{0) + 7 ±)h + 2(d + 2)h 



+ 



(d - l)a r , 
8L 2 



(d-l)I 2 + 2I 5 



= L 



(d 2 - 1) 
120 



5(0) + L 



(d-l) 
1440 



(25d 2 + 3Qd + 23) + n(lld + 5) 



(46) 



11 



+ 



L 2 \ (d-l)L 



12 



(d-i) 



+ H n + 2d 2 (H 8 + Ho 



(d-1) 



[d# 12 + 2(d + 2)H 13 + dH u ] 



L 2 
2 



(d- l)n 



12 



; (d-l) 



120 V ; 



+ L 



2 (<*-!) 

1440 



(25d 2 - 22d + 25) + 4(n + Ad - 2) 



+ L 3(^, (0)+£2 (rf-D 



120 



720 



(29d - 1) . 



(47) 



Inserting these results into Eq. (44), we obtain the partition function and all moments up to 
order e 2 oc I 2 : 



l-sL {d :} )n +e 2 L 2 



12 



(d-lfn 2 (d- l)(4n + 5d- 13)n 



288 



+ 



1440 



- <^ 3 ) 



= '-6 ,+ 



n 

72 + 



(4n + 5d- 13)n 
360(d- 1) 



Z 2 - C(/ 3 ) . 



(48) 



For n = this yields the properly normalized partition function Z = Z = 1. For n = 2,4 we 
obtain the known even moments in d dimensions 



r 



1 l9 

= 1-- + — 2 + ... , 
3 12 

„ 2/ 25d-17 l2 

= + — 7^ +••■ • 



3 90(d-l) 

In addition, we find all odd moments up to order I 2 , the lowest being 



<r> 



1 - 



I 



5d- 7 



6 180(d-l 



-/ 2 + 



1 j | 25d-4 

2 + 30(d-l) + '-- • 



(49) 
(50) 

(51) 
(52) 



Restricted to d = 3 dimensions, all moments (48) are in full agreement with the exact 
expansions of Ref. 23 obtained by solving the diffusion equation on a unit sphere. 



5 Correlation Function Up to Four Loops 

As an important test of this perturbation theory, we calculate the correlation function (7) up 
to four loops and verify its agreement with the exact expression (8). 

Starting point is the path integral representation (29) with Neumann boundary conditions 
for the two-point correlation function 

G(s,s') = (u(s)-u(s'))= [ V' d - 1 q {s)U[ q ,q';e}e^{-A^[q;e}} , (53) 

./NBC k ' 



12 



with the action of Eq. (36) for n — 0. The functional in the integrand U[q,q'} = U(q(s),q(s')) 
is the scalar product u(s) • u(s') expressed in terms of independent coordinates g' 1 as 



£/(?(*),?(*')) = u(s) -u(s') = v /l-g2 (s)v / 1 _ g2(s0+g(s)g(s0 
After rescaling — >■ y^g^ and expanding in powers of e, this becomes 

E% <?'; e] = 1 + eE^g, <?'] + e 2 f/ 2 [g, ?'] + ..., 

with 

^i[?,9l = ^i(?(s),?(s0) = ?(*)?(*0-^(*)-^V), 



^ 2 ( s )gV)-^b 2 ( s )] 2 -^[g 2 (^)] 2 



^2b,g'] = ^2(g(s),g(s')) : 

We shall attribute the integrand U[q, q'; e] to an interaction A v [q; e] defined by 

U[q,q';e] = expi K -A V [q;e]) , 



(54) 

(55) 

(56) 
(57) 

(58) 



which has to be added to the interaction action in Eq. (36) with n — 0. The small e-expansion, 
similar to Eq. (40), reads 



A v [q; e] = -eU^q, q'\ + e 2 \-U 2 [q, q'} + -U*[q, q'\ 



Thus, the perturbation expansion of the path integral (53) takes the form 



(59) 



G(s,s') 



1 



1 - ((A^[q;e] +A v [q;e})) + ^((A^[q;e} + A^[q;e})\ 



(60) 



Collecting the expansions terms of Eqs. (40) and (59) in Eq. (60) and taking into account the 
unit normalization of the partition function Z = Z = 1 reduces this to 



G(s,s') = l + £ (U 1 [q,q']) + e 2 \(U 2 [q, q'}) - (U^q, q'\ A*f[q\)d + . . . , 



(61) 



where the expectation values can be calculated, as before, using the propagator (13) with the 
Green function (33). For calculation of the rotationally invariant quantity G(s, s') = G(s — s'), 
however, the Green function A(s, s') + C must be just as good a Green function satisfying 
Neumann boundary conditions as A(s, s'). 

We illustrate this explicitly by setting C = L(a — l)/3 with an arbitrary constant a, and 
calculating the expectation values in Eq. (61) using the modified Green function. Details are 
given in Appendix C [see Eq. (152)], where we list various expressions and integrals appearing 
in the Wick contractions of the expansion (61). Using these results, we find the manifestly 
a-independent terms up to second order in e: 



(£/i[<L</]>o = (d-l)A F ( S - S ') = 



(d-1) 



s — s 



intl r 



(d-1) 
(d-1) 



-Di - -(d + l)D%- K ± -dK 2 



(62) 
(63) 



^(d-lf(s-sT. 
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This leads to the two-point correlation function 

G(s,s') = l- £ { -^±\s-s'\ +e 2 ^^(s-s') 2 -0(e 3 ) 
, \ s - s '\ (s-s') 2 

in agreement with the small-flexibility expansion of the exact expression (8) up to the second 
order in l/£. 



6 Alternative Calculation of All Moments 

It is interesting to recapitulate the result (48) by calculating all moments (6) once more in the 
way used in the previous section. To this end, we start from the path integral with Neumann 
boundary conditions 

(r n ) = J^V'^q^R^q-e] exp {-A^[q; e]} , (65) 

where the action A t( Q [q; e] is the same as in Eq. (53), while the functional R 2 [q] corresponds to 
the square of the reduced distance r 2 represented now in terms of the d — 1 coordinates q^(s) 
with the help of Eq. (54) as 

B?[q] = L- 2 [ L ds [ L ds' U{q{s), q(s')) . (66) 
Jo Jo 

Note that, in terms of the reduced distance r, the moments (65) make sense initially for all n. 
Using Eqs. (55)-(57) yields the small e-expansion 

R 2 [q; e\ = 1 + eR\[q} + e 2 R 2 2 [q} + . . . , (67) 

with 

R\\ q ] = L~ 2 ( L ds ( L ds U l {q{s),q{s')) = -L- 1 ( L dsq 2 {s), (68) 
Jo Jo Jo 

R 2 2 [q] = L- 2 [ L ds [ L ds' U 2 (q(s),q(s')) 
Jo Jo 

= -(4L)- 1 f L ds[q 2 (s)} 2 + (2L)- 2 ^ ds ^ ds' q 2 (s)q 2 (s') , (69) 
Jo Jo Jo 

where we have used the condition (23). 
The interaction associated with R n [q;e] is 

A R k,s] = -|log R 2 [q;e] = -e\R\\q\+e 2 \ [-R\[q\ + \{Rl[ q ]?} - . . . . (70) 

Adding (70) to the action A^lq; e] leads to the perturbation expansion of the path integral (65) 
up to second order in e: 



(r n ) = l-((< [g;e]+^[g;e]))o + i(Ko t [g^]+^ R [g;e]) 2 )o 



= l + e'^Rl^o + e 21 - 



(Rllq])o - (R 2 M A^[q]) + ](n - 2)((R 2 1 [q]) 2 )o 



■ (71) 
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Most of the expectation values appeared here are related through Eqs. (68) and (69) to those 
obtained before in Eqs. (62) and (63). For these, we find 



(R 2 M)o = ^ J o L ds^ds'(UMs),q(s'))) 



(d-1) 
2L 2 



[ L ds f 1 
Jo Jo 



ds' I s — s' 



-L 



(d-1) 



(72) 



(R 2 Mo-(Rl[q}^[q}) = ± /^ Sj ( L rf S '[(f/ 2 (g( S ),g( S ')))o-(f/i(^),g( S / ))<o 1 M)o 

( ( L ds >( s - S >)l = L 2(^-l) 2 

Jo 



48 



(73) 



The last expectation value in Eq. (71) is calculated as 



' fL 'ds [ L ds'(q 2 (s)q 2 (s')) 
Jo 



L 2 Jo Jo 
(d-1) 



L 2 



(d - i)ii + 2/ 5 



= L 



2 (<Z- l)(5d- 1) 



180 



(74) 



Inserting the results (72)-(74) into Eq. (71), we obtain the same perturbation expansion for all 
moments up to order I 2 as in Eq. (48), but in the slightly different form 



<r»> = l-eL {d :} )n + e 2 L 2 



12 



(d-l) 2 n (d-l)(5d-l)n(n-2) 



96 



+ 



1440 



0(s 3 ) 



n (5d — l)n(n — 2) 
24 + 360(d- 1) 



I 2 - 0(l 3 ) . 



(75) 



7 Radial Distribution up to 4 Loops 

We now turn to calculation of the most important quantity characterizing a polymer, the 
distribution function (25). Going over to the Fourier transform of the 5-function, which en- 
forces the reduced end-to-end distance r = L^ 1 J L ds \jl — q 2 (s) in Eq. (25), yields the Fourier 
decomposition 



P(r; L) = S? J £ e-^-D P ( k; L ) , 



(76) 



where P(k; L) has to be calculated from the path integral with Neumann boundary conditions 

P(k;L) = J^ c V' d - 1 q (s)exp{-A^{q;e}} . (77) 
The resulting total action 

A%[q] = A[q] + A™[q bl q a ] - (ik/L) £ ds (y/l - q 2 (s) - l) , (78) 
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after rescaling the coordinates q^(s) — > \/eq tJ '(s), reads explicitly 



* 2 



1 — eq 2 



+ \ 6(0) log(l - a? 2 ) - ^ (yfl - eq- - l) I 



' [ q 2 (0) + q 2 (L)]-eL^ A°[q}+A^[q;e}. 



(79) 



As before in Eq. (40), we expand the interaction in powers of a small coupling constant e. 
The first term coincides with Eq. (41), except that p n {s) is now replaced by Pk(s) = 5k + 
[S(s) + 5(s - L)] /2 with 5 k = 5(0) - ik/L: 



ds - 
2 



(q(s)q(s)) 2 -p k (s)q 2 (s) 



L— . 



(80) 



The second expansion term «4™ t2 [(z] is simpler than that in Eq. (42) by not containing the last 
nonlocal contribution 



s )?-\(m- l ^\ q\s)\q%s). 



(81) 



Apart from that, the perturbation expansion of the path integral (77) has the same general 
form as in Eq. (44) and, therefore, reads 



P(k;L) 



(A m k t [q;e]) + -(A 1 *[q;e} 2 ) 



1 - e (A^[q]) + e 2 -« t fc 2 [g])o + -AA^[qf)o - 



(82) 



The expectation values can be expressed in terms of the same integrals listed in Appendices A 
and B as follows: 



<4r[?]>o = 



(d-i) 



*h + dI 2 -\±(W)-\±(L,L) 



L— 

8 



—L 



(d-1) [(d-l)-ik] 



12 



(83) 



(d 2 - 1) 



ik 



5(0) + —jh + 2(d + 2)I 4 
L - —5(0) + L 2 



4 

(d 2 -l 



(d 2 - 1) [7(d + 2) + 3ik] 



120 



720 



(84) 



12 



d 



m + 2L - {h + T 



R 
~8 



+ 



(d-1) 



+ H n + 2d 2 (H s + H 9 ) 



if* - 2(H k 2 + Hi - H 5 ) + H 6 - Ad(H k 4 - H 7 - H w ) 
(d-1) 



+ 



L2 (d-l) 2 [(d-l)-tk] 2 
2 ■ 12 2 



120 



1440 



+ L ^ d :l )d 5(0) + L 2{d - l) 



120 



720 



[dH 12 + 2(d + 2)H 13 + dH 14 ] 

(13d 2 -6d + 21) + Aik(2d + ik) 
(29d- 1) . 



(85) 
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In this way, we find the large-stiffness expansion of the path integral (77) up to order e 2 : 
P(k;L) = 1 + eL 



(d_1) [{d-l)-ik]+s*I? {d - l) 



12 



x 



1440 

+ (d - l)(5d 2 - lid + U)} + 0(e 3 ) . (86) 



{ikfihd - 1) - 2ik{5d 2 - lid + 8' 
The Fourier transform (76) of the expansion (86) yields the end-to-end distribution function 

P(r;l) = S?{6(r - 1) + l - [5>(r - 1) + (d - 1) 5(r - 1)] + ^ _ ^ [(5d - 1) 5"(r - 1) 

+ 2{5d 2 -lld + 8)5'{r-l) + (d-l){5d 2 -lld+U)5{r-l)]+0{l 3 )}. (87) 
This representation is convenient to calculate the moments 

( r ») = S d [ drr^-V P(r;l). 



(88) 



Indeed, inserting the distribution function (87) into Eq. (88) yields directly the moments (48) 
found before by the independent calculation. 

To make the contact with the results of Ref. 23 derived in d — 3 dimensions, we rewrite the 
distribution function (87) in the form 



P(r; /) = — \p (l) S(r - 1) + Pl (l) 5'(r - 1) + p 2 (l) 5"{r - 1) + p 3 (l) 6"'(r - I) + . . .] 

where the coefficients Pi(l) are expanded for d = 3 up to order I 3 

3t 



Po(l) 
Pi(0 

P2(0 

Ps(0 



l + i/ + i^/ 2 + -- 

3 180 2 3 -315 



,o 2 X3 o 38 o 

' + --= 1 + 3' + 45' +315* + -' 



/ 3 + ... = - t+-t 2 + 

3 9 315 



1,1 ,2, 34 „ . 1 . 2 . 34 

6 + 18 + 23-315 

7 1 7 1 

/ 2 + — / 3 + ... = — t 2 + — t 3 + . . . 

360 2 3 -21 90 21 



t 3 + ... , 



31 



2 3 • 1890 



Z 3 + ...= 



31 

1890 



t 3 + . . . . 



(89) 

(90) 
(91) 
(92) 
(93) 



Inserting I — 2t in the last lines of Eqs. (90)-(93) shows the agreement with the end-to-end 
distribution function G(r; t) of Ref. 23 reduced to the form (89) after re-expansion in Appendix D. 

Let us make now the contact with the distribution functions of Ref.. 25 To this end, the 
expansion (86) can be reassembled as 



P(k;L)=P (k;L) U + ^-Jll + 



(d-3) (5d 2 -lld+14) 



180(d- 1) 



360 



I 2 + 0{l 3 ) , (94) 



where we have factored out the series containing the expansion terms with the same powers of 
e, k and L, 



P {k; L) = 1 - ^— -iikeL) + -A^ '-{iksVf - 



2 2 -3 



2 5 • 3 2 • 5 



(95) 
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By setting Cj 2 = ikeL, we identify this with the beginning of the expansion of a functional 
determinant 

(d-l)/2 



P (w) = 



a; 



sinhu) 



= 1 



d-l.o d-l(5d-lk 4 

cu + -u; 4 + . . . 

12 1440 



(96) 



valid for Co small. One can replace, however, the expansion (95) by the exact expression (96) 
for P (u>) and integrate over Cj in the Fourier transform (76). This would ignore the smallness 
of Cj. In such a way, the distribution function with the only determinant (96) was calculated 
for d = 2, 3 in Ref. 25 from the u)-integral 



P (r;l) oc 



i r°° duj 



3 -[i 2 (d-l)(r-l)/2Z 



2ni 



CjP (Cj) , 



(97) 



which was found to give, after a proper normalization, good radial end-to-end distribution func- 
tions for large stiffness. It represents the leading contribution to the distribution function (76), 
since Pq{Cj) enters Eq. (94) as a prefactor with no dependence on the flexibility I. To find 
corrections, the rest of expansion (94) in powers of small / must be known for all Cj. This can 
be done by reorganizing the perturbation theory as follows. 

The determinant (96) can directly be obtained from the path integral of a simple harmonic 
oscillator with Neumann boundary conditions 15 : 



Zl = [ V' d - l q{s)e' A °^ = P {Cj), 

JNBC 

where -A ^] is the harmonic action 



ds 



q 2 (s)+cj 2 q 2 (s) 



(98) 



(99) 



With the redefinition uo 2 = ike/L, it represents a "free" part of the total action (79), whose 
expansion in powers of small e becomes 



while the first expansion term reads 

"I 



ds 



(qq) 2 -l5(0)q 2 +lu 2 (q 2 ) 2 



1 

4 



q 2 (0) + q 2 (L) 



R 



(100) 



(101) 



Note that the interaction (101) contains the harmonic terms even at lowest order in e. 

The result (98) represents the leading term in the large-stiffness expansion for the path 
integral (77). By expanding the path integral (77) with respect to the harmonic integral (98), 
we find the higher-order fluctuation corrections. This yields 

P&L) = Zl (l-eMSfWo + .-O , (102) 

where the expectation values are evaluated using the basic propagator (13), which contains now 
a harmonic Green function of the unperturbed action (99). With the zero-mode subtracted, 
this Green function explicitly reads 



, 2 ^ cos(nns/L) cos(mrs' /L) 



n=l 



u 2 + n 2 n 2 /L 2 



coshcj(L— | s — s' |) + coshcj(L — (s + s')) 1 



2u sinh ujL 



Luj 2 ' 



(103) 
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As a result, we obtain the next-to-leading contribution 



8 4 \ Cj Cj 2 



r id 2 — 1) / , 9 , cothcj 
+ L - - coth 2 £ + 



32 



2cD 2sinh Cj 



(104) 



Substituting now the leading (98) and the next-to-leading (104) terms into expansion (102) 
yields 



P(Cj; 1) = P (Cj) R(cu; I) , 
where the factor R{Cj; I) for all Cj has the small-/ expansion 



(105) 



R(cj; 1) = 1 + 1 



(3d -5) Ad+1) (cothu 1 \ 1 /cothw 1 



16 



32 y Cj sinh a)/ 2 y Cj Cj 



^ +o(n- ace) 



Expanding (105) in powers of uj for Cj small yields, of course, the expansion (94) to the order 
of I. Nevertheless the analytic form (105) is convenient for the use in the Fourier representa- 
tion (76) in terms of Cj: 



P(r; I) = S d l r ^ e -^(--i)(r-i)/2i QP{Qj . z) 



(107) 



An important observation for performing the u)-integral (107) is now that the integrand (105) 
can be expressed in terms of Pq(oj), and the derivatives Pq(Cj) and Pq{Cj) with respect to Cj. To 
show this, we set 



Po(Cj) 



cj 



(d-l)/2 



sinh u> 



,-/(£ 



with 



/(a)) = (log sinh Cj — logcD) . 



In terms of /(u>), the factor (106) becomes 

R(r,-n i i (3d-W (d-15)l f'(u) (d+l)l „ 

Substituting Eqs. (108) and (110) into Eq. (105), we use the identities 
-f'(u)e-fW = P>(Cj), 



(108) 



(109) 



(110) 



-f\0j)e 



-m 



2(d-l) 



(d-iy 



This brings, finally, the expansion (105) in the form 

\ , (5d-9)f 



Po(^) 



(111) 



P{Cj-1) 



32 



(d+13)Z Z 
P oM - 7777 — 77^ P oM + 



16(d- l)w 



8(d- 1) 
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With the representation (112), the integration of Eq. (107) is straightforward. We substitute 
the binomial series 



(d-l)/2 



sinhu) 



(2a>) 



(d-l)/2 



E( 

fc=0 



-(d-l)^ 



-(2fc+(d-l)/2)w 



(113) 



and express the resulting integrals in terms of a parabolic cylinder functions. In the physical 
case of three dimensions, where a) 2 = ikl, the answer reads 



P(r;l) 



E 



^V^lti [(l-r)//] 3/2 



cxp 



n 



l/2f 



(l-r)/l 



x 



1 + 



n+ 1)Z 



n- 1/2 

;i-o/z. 



+ 3 ( n-l/2)^_ (i _ r) 
v2 



(114) 



where the second Hermite polynomial is H%(x) = 4x 2 — 2. The first term in the brackets 
represents the leading contribution found before in Ref., 25 the other are the new next-to-leading 
corrections. The proper normalization is ensured automatically in each perturbative order / 
in the small-flexibility expansion. Since the next-to-leading terms are suppressed as I ~ 1/k, 
there is no drastic change in the behavior of the radial distribution with respect to the leading 
contribution. 



8 Conclusion 

In conclusion, we have developed a systematic perturbation theory capable of yielding the 
polymer properties near the rod limit from a path integral formulation. This has yielded the 
large-stiffness expansions for the experimentally relevant end-do-end distribution, all the even 
and odd moments, and the correlation function in d dimensions as power series in the flexibility 
/ up to the order I 2 . All these results have been obtained without the use of the diffusion 
equation on a unit sphere. In subsequent work we shall use variational perturbation theory 15 ' 41 
to extend the results from small to large flexibility, i.e., to find results connecting the above 
results smoothly with the random-chain limit. 



Appendices 

A Basic Integrals 

The following well-defined integrals appear throughout the calculation of the Feynman dia- 
grams: 

A(0,0) = A(L,L) = |, (115) 

h = [ L dsA(s,s) = ^-, (116) 
Jo o 
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h 
h 
h 
h 

A(0,L) 

h 

h 

h 

h 

ho 
'A (0,0) 

hi 
In 



fL 

/ ds A(s, s) A 2 (s, s) 
Jo 

/ L rfs [ L ds' A 2 (s,s') -- 
Jo Jo 

A(L,0) = -|, 



7L 2 
" 360 

90 ' 



A 2 ( s ,0) + A 2 ( s ,L) 



2L 3 

"45"' 
L 3 



JO 

= ( L ds T ds'A(s,s)A\s,s') = ±-, 
Jo Jo 45 

= / L [ L ds' A (s, s) A(s, s') A (s, s') 

JO JO 

= f L dsA(s,s) pA 2 (s,0) + A 2 (s,L)" 
Jo L J 

= [ L ds A (s, s) [A(s, 0) A (s, 0) + A(s, L) A (s, L)] 
Jo 



180 ' 
ILL 2 

90 



17L 2 
360 ' 



■A(L,L) = --, 



[ L ds [ L ds' A(s, s) A (s, s') A (s', s') 
Jo Jo 



360 ' 



f ds r ds'A(s,s')A 2 (s,s') = ^-. 
Jo Jo 90 

B Loop Integrals 



(117) 
(118) 
(119) 
(120) 
(121) 
(122) 
(123) 
(124) 
(125) 
(126) 
(127) 
(128) 
(129) 



We list here the Feynman integrals evaluated with dimensional regularization rules whenever 
necessary. In the calculation, they occur either from the expectations (45)-(47), or from the 
expectations (83)-(85). Accordingly, we encounter the integrals depending either on p n (s) = 
S n + [S(s) + 5(s - L)] /2 with 5 n = 5(0) - a n /L, or on p k (s) = 5 k + [5(s) + 5(s - L)) /2 with 
S k = 5(0) - ik/L: 



O H[ l{k) 


= J Q ds J ds'p n(k) (s)p n{k) (s')A 2 (s,s') 

= 5 2 n{k) h + S n{k) h + \ [A 2 (0, 0) + A 2 (0, L) 






(130) 


oo H 2 k) 


= ( L ds [ L ds' Pn(k) (s')A(s,s)A 2 (s,s') =5 
Jo Jo 


n(k) h + y, 




(131) 


CO H? k) 


= J ds J ds'p n(k) (s'YA(s,s)A 2 (s,s') = 


5 n (k) h + y 


5(0), 


(132) 


OO Hl {k) 


= £ dsJ L ds'p n(k) (s')A (s, s)A (s, s')A(s 


s') = 5 n ( k ) It 


.ho 


(133) 
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When calculating Z n , we need to insert here 5 n = 5(0) — cr n /L, thus obtaining 

= ^5 2 (0) + ^(3-d-n)5(0) + ^[(A5-2Ad + Ad 2 )-An(6-2d-n)],(13A) 

HI = ^5(0) + ^(15-Ad-An), (135) 

1?J = ^5 2 (0) + ^(3-d-n)5(0), (136) 
L 3 L 2 

H? = 5(0) + (21 -Ad -An), (137) 

4 180 V ; 720 V ; ' V ; 

where the values for n = correspond to the partition function Z = Z . The substitution 
5k = 5(0) — ik/L, required for the calculation of P(k; L), yields 

7-4 r3 t2 rrj2 

Hi = ^S 2 (0) + -[2 + (-ik)]5(0) + -(-ik)[A + (-ik)} + —, (138) 

H k 2 = ^*(0) + ^[ll + 4(-ifc)], (139) 
H* = ^5 2 (0) + ^[2 + (^k)]5(0), (140) 

H\ = 7^°) + ^[ 17 + 4 H^- (141) 

The other loop integrals are 

rL rL ^ 



GGOH 6 = J dsj^ ds'A(s, s yA 2 ( s ,s'yA(s',s') = 6(0)1? = ^5(0), (142) 

COD H 6 = £ ds£ ds"X(s, s)A 2 (s, s'yA(s', s') = S 2 (0)I 5 = |j <5 2 (0), (143) 

000 H 7 = £ 'da £ d S 'A( S , S yA( S ,s')A(s,s , yA(s , ,s , ) = 6(0)h = ^5(0), (144) 

000 #8 = £ ds£ ds'A (s, s)A (s, s') A(s, s')A (s', s') = (145) 

(146) 

O'X)h 10 = £ d s £ d s 'A( s , s yA( s , s yA( s , s yA( s ',s') = |- ^ = (w?) 

Q00#n = [ L ds t ds'A(s,syA 2 (s,s')A(s',s') 
Jo Jo 



H 12 = ( L ds [ L d s / A 2 (s,s')A 2 (s,s') = ^-, (14!)) 
Jo Jo yu 

<E> H 13 = ( L ds ( L d8'A(s,s , yA(8,s')A(s,s'yA(s',s)= 1 -±-^-^ = -^-, (150) 



ho ^8 

2 L 


-H s = 


7L 2 
360 


5 


/9 h 


7L 2 






4 2L 


" 360 






h|) 2 - 


2(/ 3 - 
L 




- 2/ 4 


= s>> 


L 2 






h hi 






L 2 


2 2L 


2~~ 




720' 



H 14 = ds ds , A 2 (s,s , yA 2 (s,s') =6(0)I 3 - 
Jo Jo 



L lL 2(I 3 -I 12 ) h 

2/ 4 + - 

LK,„ IIL 2 



L 2h + f^ 

lQ Jo L, L, 



(148) 



+ 2 [A 2 (L,L) A(L,L) - A 2 (0, 0) A(0, 0)] - 2H 13 = ^6(0) + (151) 
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Calculating the integrals (146)— (151) required the regularization rules of Refs.. 37 ' 38 To compute 
these unambiguously, we must first use partial integrations together with Neumann boundary 
conditions so that we can apply subsequently the equation 'A (s, s') = A"(s, s') = l/L — 5(s, s'). 
In this way, most of the integrals can be expressed in terms of the basic integrals in A. For the 
integrals (148) and (151), the above procedure had to be applied twice. 



C Integrals Involving Green Functions with Arbitrary 
Constant 



To demonstrate the translational invariance of results showed in the main text we use the 
modified Green function 



[S,S) 3 2 2 + 2L 



(152) 



with an arbitrary constant a. The various relations fulfilled by this Green function are listed 
below assuming s > s', for brevity. 



The useful relations are 

A F (s, s') = A F (s - s') = A(s, s') - ±A(a, s) - ±A(s f , s') 



-2 [S ~ S 



L>i(s,s') = A 2 (s,s') - A(s,s)A(s',s') = (s - s') 



s(L-s) (s - s')(s + s') 2 La 



+ 



4L 2 



D 2 (s,s') = A(s,s)-A(s',s>) 
The useful integrals are 



{s - s')(s + s' - L) 



(153) 
, (154) 

(155) 



JJ s ,s')= / dtA(t,t)A(t,s)A(t,s') 
Jo 



(s 4 + s' 4 ) (2s 3 + s' 3 ) 



4L 2 



3L 



((a + 3)s 2 + as' 2 ) sa (20a - 9) 2 

H — — Lj -h Lj , 

6 3 180 



w n f L , -a / \ a / \ • a / a (-2s 4 + Qs 2 s' 2 + 3s' 4 ) 
J 2 (s, s') = dt A (t,t)A(t, s) A (t, s') = ± 



+ 



o 24L 2 

(3s 3 - 3s 2 s' - 6ss' 2 - s' 3 ) (5s 2 - 12ss' - 4as' 2 ) s'a (20a - 3) 



12L 



J 3 (s,s')= [ L dtA(t,s)A(t,s') 
Jo 



24 -T L + 

[s 4 + 6s 2 s /2 + s' 4 ) (s 2 + 3s' 2 )s 



720 



60L 



+ 



(s 2 + s' 2 ) (5a 2 - 10a + 6) 
-L H 



L 3 . 



6 45 
These are the building blocks for other relations 



(156) 



L\ (157) 



(158) 



(U2(s,s')) = 



(d-i; 



Dl (s,s')-^±^D 2 2 (s,s') 



(d-l)(s-s') 



2La 



(d-3)s-(d+l)s' _ (d-l)s 2 - (d+l)s' 2 d(s - s')(s + s') 2 
+ 2 L + 2L 2 



(159) 
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K x {s, s') 
K 2 (s, s') 

K 3 (s,s') 

K 4 (s,s') 
K 5 (s,s') 



Ji(s, s') - ^Ji(s, s) - ^Ji(s', s' 



— (s — s 



La (s + s') (s 2 + ss' + s' 2 ) 
~6 4 + 6L 



(160) 



J 2 (s, s') + J 2 (s', s) - J 2 (s, s) - J 2 (s', s 
-(*- s') 2 



(s + 2s') (s + s'f 



1 (5s + 7s') (s + s') 2 
4 12L + 4L 2 

1 i 

J s(s, s') - - J 3 (s, s) - - J 3 (s, s') = -(s - s') 2 

A(0, s) A(0, - iA 2 (0, S ) - iA 2 (0, = - (S " S/)2(S 8 + S/ ' 2L)2 , (163) 



6 



8L 



(161) 



,(162) 



A(L, s)A(L, sf) - \A 2 (L, s) - U 2 ^ s>) = 8 f ? ■ 



(164) 



D Relation with Distribution from Diffusion Equation 



The distribution function derived in Ref. 23 for d = 3 dimensions has the form 



fo(t) S(r - 1) + Ml 5\r - 1) + ^ 5"{r - 1) + ^ 5"'(r - 1) + 



,(165) 



where the coefficients fi(t) are the polynomials in powers of t. To order t 3 , their expansions 
read 



/„(*) = l + it + lt 2 + A,3 + ._ 
hit) = hit) -I, 



90 630 



/s(*) 



31 

1890 



■t a - 



(166) 
(167) 
(168) 

(169) 



Depending on the reduced distance r, the expression (165) cannot be directly comparable with 
the distribution function (89). To remove this, we expand 



= 1 - (r - 1) + (r - l) 2 - (r - l) 3 + . . . , 



r (r - 1) + 1 

Substituting the expansion (170) into Eq. (165), we use the properties of 5-function 



(170) 



x5'(x) = -5(x),x n 5'(x) = 0,n>l, (171) 
x5"(x) = -25'(x),x 2 5"(x) = 25(x),x n 5'(x) = 0,n>2, (172) 
x5"'(x) = -35" (x) , x 2 5"'(x) = 65' (x) , x 3 5'"(x) = -65(x) , x n 5'(x) = ,n > 3 . (173) 

This brings the distribution function (165) in the form (89) with 



Mt) = 2 /o (t) + 2^ + 6^ 



2 13 o 38 q . _ . 

1 + ... = 1 + _ t + _ (2 + _ t 3 + ..., 
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Pl (t) = /o(i) + 2^+6^-l + ... = it + |t= + |±t 3 + ... , (175) 

m(t) = ^r + 3 ; iP + '-- = s i2 + 2T t3 + "' ' (176) 
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